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Abstract 

We propose a quantum approach to nonequilibrium dynamics which 
combines the successful aspects of classical-statistical simulations on 
a lattice with the ability to take into account quantum corrections. 
It is based on the 2PI effective action for inhomogeneous fields and a 
volume average. This procedure does not involve any double count- 
ing which could appear in sampling prescriptions for inhomogeneous 
quantum evolutions. As an example, we study nonequilibrium dy- 
namics of defect formation in 1 + 1 dimensional relativistic scalar field 
theory and compare to insufficient descriptions based on homogeneous 
quantum fields. The latter are analyzed in detail by coupling the field 
to an external source, such that the emerging infiuence of defects can 
be studied by lowering the source to zero. 



1 Introduction and outline 



Defect formation in symmetry breaking phase transitions is a nonequilibrium 
process with important apphcations for condensed matter physics or even 
early universe dynamics [1]. It involves nonperturbative djTiamics far from 
equilibrium which is a challenge for theoretical descriptions of quantum many 
body systems. Often the physics is approximated by classical-statistical field 
theory, which can be simulated on a lattice using Monte Carlo methods 
and numerical integration [2]. Classical Rayleigh- Jeans divergences and the 
lack of genuine quantum effects, such as the possibility to approach thermal 
equilibrium characterized by Bose-Einstein or Fermi-Dirac distributions, 
limit their use [3]. Lattice simulation methods for nonequilibrium quantum 
field theories, such as based on stochastic quantization [1], are still in its 
infancies and it would be very desirable to be able to compute the impact of 
quantum corrections on defect formation using alternative approaches. 

There has been substantial progress in recent years in the quantitative 
understanding of the nonequilibrium time evolution of quantum fields using 
approximate functional integral techniques. All information about a given 
nonequilibrium quantum field theory can be encoded in an effective action, 
which is the generating functional for all correlation functions. While 
naive perturbative expansions lead to secular behavior and fail to describe 
the time evolution, suitable approximations may be efficiently based on 
the two-particle irreducible (2PI) or higher nPI effective actions [3]. The 
nonperturbative 2PI 1/N expansion to next-to-leading order (NLO) [5], 
where denotes the number of field components, has been particularly 
fruitful. This approach was applied to a variety of problems, related to 
thermalization El E], inflationary preheating [S], fermion dynamics [H], 
transport coefficients [10], nonthermal fixed points [11] and critical exponents 
of phase transitions [12], or cold atoms [13]. The 2PI 1/N expansion has 
also been implemented in classical-statistical field theories, where it can be 
tested against full simulations on a lattice. Remarkable agreement was found 
between simulation results and those from the NLO approximation in various 
dimensions for not too small N [2]. Also NNLO corrections have been 
addressed [T5] . 

In contrast to those successful applications, it was pointed out in Ref. [16] 
that the 2PI 1/A^ expansion to NLO does not reproduce known results of 
phase transition dynamics for A^ = 1 and N = 2. Of course, a 1/A^ expansion 
is not expected to work well for too small A^ and quantitative deviations were 
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Figure 1: The dynamics of the spatially homogeneous field expectation value 
(f)(t) in the presence of defects can be replaced by a description using an 
inhomogeneous field x) and a volume average. This allows one to study 
topological defect formation in quantum theories using inhomogeneous 2PI 
1/N or loop expansions. 



observed before in this case [H]. However, the results of Ref. [16] suggest that 
the discrepancies arise because topological defects are present in this case, 
which are not properly described. This suggestion follows a long-standing 
discussion about the question of whether approximate functional integral 
techniques can describe topological defects. 

In this work we follow the nonequilibrium time evolution of the order 
parameter (f){t) and two-point correlation functions in a scalar 0{N)- 
symmetric field theory with quartic self-interaction. This extends the 
study of Ref. [12], where defect formation from the dynamics of two-point 
correlation functions with = is considered. We further couple the field to 
an external source J. For large enough J topological defects are suppressed 
and we can study their emerging infiuence in detail by lowering J for given 
N. We find that for not too small source there are only relatively small 
deviations between classical-statistical simulation results and those obtained 
from the 2PI expansion to NLO even for small A^. On the other 

hand, as J becomes sufficiently small such that topological defects become 
relevant, we find significant deviations. Our results confirm the conclusion 
that contributions from topological defects present for N = 1 and N = 2 are 
not properly described at NLO in the 2PI 1/N expansion. In general, we 
find that NLO results for N = 2 agree somewhat better to those obtained 
from classical-statistical simulations than for = 1 as expected for a large- 

expansion. It would be very interesting to extend the study to NNLO 
though the computational effort becomes considerable [15] and alternative 
approaches can be more appropriate. 
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Here we propose an alternative quantum treatment, which allows one to 
quantitatively describe topological defect formation already at low orders 
in a 1/N expansion. At leading order the approach is similar to the 
inhomogeneous Hartree approximation considered in Ref. [l7j. It constructs 
the spatially homogeneous quantum system from a corresponding theory 
for an inhomogeneous "foreground" field </'(t, x) and a volume average as 
indicated in Fig. [T] The additional spatial x-dependence of the field x) 
is constrained by the fact that its volume average corresponds to the spatially 
homogeneous field expectation value of the original quantum theory: 

0(t) = (0(t,x))y, (1) 

where {■ ■ ■)v denotes an average over some local volume V whose size is 
determined by the statistical error one wants to achieve. Similarly, the 
higher n-point functions of the inhomogeneous description are related to 
those of the homogeneous quantum system by a volume average. We 
emphasize that without approximations the volume averaged inhomogeneous 
description agrees to the original system. In particular, this procedure does 
not involve any double counting which could appear in sampling prescriptions 
for inhomogeneous quantum evolutions. Here the description only exploits 
the fact that a spatially homogeneous system is by definition translation 
invariant in space. We show that the classical-statistical approximation is 
included already in the inhomogeneous quantum approach to leading order 
in the 2PI expansion, while quantum corrections enter at next-to- leading 
order. 

The outline of the paper is as follows. In Sec. [2] we describe classical 
kink and texture solutions in one spatial dimension for an 0{N) symmetric 
real scalar field theory. In Sec. [3] we include classical-statistical fluctuations 
by numerically integrating the classical field equations and Monte Carlo 
sampling. Sec. H] introduces the quantum field theoretical description using 
the 2PI effective action. We discuss the role of quantum corrections and 
describe the classical-statistical limit. We consider in Sec. EJthe conventional 
approach, where the 2PI effective action is evaluated for a time-dependent but 
spatially homogeneous field expectation value 4)(t). The results are compared 
with those from classical-statistical simulations for = 1 and N = 2 for 
different J. In Sec. Owe introduce our inhomogeneous approach and compare 
with classical-statistical simulations. Finally, in Sec. [7] we conclude and give 
an outlook. 
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2 Topological defects in the 0{N) model 



We consider an 0{N) symmetric A^-component real scalar field theory with 
quartic interaction in 1+1 dimensions following Ref. [16]. In addition, we 
couple the classical fields ipa(t,x) (a = 1, . . . , with t denoting time and 
X is the variable for the spatial dimension) to explicitly symmetry breaking 
source terms Jq. The classical field equation of motion reads 

d^ifa + idtifa - dl^Pa + /^Va + ^j{.VhVh)Va - Ja = , (2) 

oiV 

with mass parameter /i^ and quartic coupling A. Summation over repeated 
indices is implied. The field equation of motion ([2]) contains a damping term 
•ydfipaityx) whose relevance will be discussed below. 

The non-derivative part of ([2]) can be obtained from a quartic classical 
potential with linear source term, 

Vi'P) = ^f^^^a^a + -^{^a^af - Ja^a ■ (3) 

For the study of topological defect formation the initial classical potential is 
taken to be of free field form, 

initial (V^) = ^m^<-Pa<^a, (4) 

where > 0. At time t = the potential changes instantaneously to 
03]) with /i^ = — < 0, so the system is "quenched" into the phase with 
spontaneous symmetry breaking. For A^ = 1 the classical potential then 
exhibits two degenerate minima at ip = ±v with 




in the phase with spontaneous symmetry breaking for J = 0. This is shown 
graphically in Fig. |5J The topological defect or kink solution interpolates 
between these separated minima. The stationary classical solution of ([2]) 
reads 

X 

V5kmk(a;) = wtanh- , (6) 

where d = \/2/\fi\ is the kink thickness. Below, we will study the behavior of 
kink solutions also for no n- vanishing source J. From Fig. |2]one observes that 
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Figure 2: Classical potential as a function of the field for different values 
of the source J for = 1. For later comparison with simulations in Sec. El 
in lattice units of we use /i^ = —0.49 and A = 0.6. Rising J from zero lifts 
the degeneracy until only one minimum remains for J > J* ^ 0.42. 



by sufficiently increasing J the potential starts to exhibit only one minimum. 
This happens around the critical value J = J* with 




Since kink solutions interpolate between two minima of the potential, they 
will be absent for J ^ J*, which is discussed in more detail in Sec. [3] on 
classical-statistical simulations. 

For N = 2 the minima of the classical potential in the phase with 
spontaneous symmetry breaking lie on a circle in the absence of sources. 
The classical equations of motion have texture solutions, which wind around 
the vacuum manifold and become indistinguishable from the vacuum in the 
infinite volume limit. Writing the two real field components as the sum of a 
real and imaginary part, texture solutions are 

^r{x) + i^i{x) = V e2-*^-^/^ , (8) 

where the winding number N^j is an integer. One observes the dependence 
of the texture solutions on the volume L. For N > 2 there are no topological 
defects for our model. 
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Finally, we will be interested in observables which are computed from 
ensemble averages. The classical-statistical ensemble average is obtained 
from a normalized initial distribution P[ip^,TT^] with </3°(x) = ipa{t = 0,x) 
and 7r°(x) = dt(pa(t,x)\t=o such that 

{Mt,x)),i = y"P(/.WP[/,7r>,(t,x) (9) 

for the one-point function and accordingly for any ra-point function. The 
measure indicates integration over classical phase space, J Vip^Vir^ = 
I Y[aY[x^'^ai^)d''^ai^)- cousidcr spatially homogeneous ensembles 

relevant for early universe inflaton dynamics, such that the field average 

(paA't) = {<^a(t,x))c\ (10) 

depends on time only. Similarly, the connected two-point function 

Fab,cl{^^ X-y) = {'fa{t, X)(fb{t', y))cl - 0a,cl(i^)0f,,cl(^') (H) 

depends only on the relative coordinate. In Sec. [3] we will construct these 
ensemble averages from many individual solutions of ([2]) for initial conditions, 
which are sampled from a Gaussian initial distribution. We also note that 
the damping term '-fdt^p in the equation of motion ([2]) is needed in order to 
be able to observe a non- vanishing field average flTUl) at late times, because of 
the absence of spontaneous symmetry breaking in a 1-1-1 dimensional system 
above zero temperature. 

Since we will be interested in correlation functions, it is important to 
know how topological effects manifest themselves in correlators. A qualitative 
understanding may be obtained from considering the special case of randomly 
distributed kinks, which for = 1 may be represented in Fourier space 

by m 

Here n denotes the number density of kinks, which will be used as a fit 
parameter in analyzing numerical results from actual classical-statistical 
simulations. Correspondingly, for = 2 we employ [16] 

FtcKtuvcip) = V ^ —^exp \-^p 1, (13) 

where v is taken to denote the same value as in (l5j) and the length scale ^ 
will be taken ClS db fit parameter. 
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3 Classical-statistical simulation 



In this section we consider the classical field theory on a spatial lattice. In 
the following, all quantities are expressed in appropriate units of the spatial 
lattice spacing a^. If not stated otherwise we choose a lattice with length 
L = 256. For the time evolution we discretize the time in steps of typically 
At = 0.1. The other relevant parameters are the mass parameter /i^ = —0.49, 
coupling constant A = 0.6 and damping rate 7 = 0.4. 

We sample the initial conditions for the classical equation of motion (|2]), 
where the initial fields ipa(t = 0,p) = Aa{p) exp{iaa{p)) are taken in spatial 
Fourier space from a Gaussian ensemble for Aa{p) with a random complex 
phase aa{p) respecting ipa{t,p) = ipl(t,-p), and similarly for dtipa{t,p)\t=o- 
As a consequence, the initial field average (fTOj) and also its initial time 
derivative is taken to vanish: 

0a,cl(t = O)=O , dt(PaAt)\t=0 = . (14) 

The initial two-point function f lTT]) and derivatives are of free-field form. In 
Fourier space they read 

Fab,cl(t,t';p)\t=t'=0 = Sab, 

Z(jJp 

dtFabAt,t']p)\t=t'=o = 0, (15) 

dtdt'Fab,cl{'t^t''^P)\t=t'=0 = Sab , 

where oOp = ^/p'^ + frfi with = — yU^. We observe good convergence by 
sampling the initial conditions over 1000 runs and exploit spatial translation 
invariance to further improve statistics by taking lattice averages. 

We first compare the classical-statistical evolution for the undamped 
(7 = 0) and damped (7 > 0) case for = 1. Fig. [3] shows the classical 
field average 0ci given by (fTOl) as a function of the homogeneous source J at 
late times. As expected, in the absence of damping no spontaneous symmetry 
breaking with (pci 7^ is observed in the limit J — t- 0. In contrast, the result 
with damping (7 = 0.4) exhibits spontaneous symmetry breaking and defect 
formation may be observed in this case. We find that for the considered 
parameters the field expectation value is rather accurately described by the 
value of the minimum of the classical potential v ~ 2.21 given by ([5]). 

We are particularly interested in the classical-statistical correlation 
function flTT]) at equal time, t = t' . According to ffT2]) . a stationary value 
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Figure 3: The classical field average 0ci as a function of the source J at late 
times for = 1. As J goes to zero a non- vanishing (p^i is observed for the 
damped case (7 = 0.4), while (f)c\ approaches zero in the absence of damping 
(7 = 0) as expected. For comparison, the solid curve shows the solution 
obtained from minimizing the classical potential (|3]). 



for this quantity in Fourier space should give direct information about the 
kink density. The upper graph of Fig. H] shows for = 1 the spatially 
Fourier transformed correlation function for low momentum, p ^ 0, as a 
function of time for different values of the source J and 7 = 0.4. After an 
initial fast growth, one observes a comparably slow evolution for small J. 
With increasing source this quasi-stationary period is diminished and clearly 
absent for J > 0.4. The latter value is close to J* given by ([7]) for which the 
classical potential no longer exhibits two distinct minima (see Fig. [2]) and 
defect formation is not expected in this case. 

The lower graph of Fig. H] shows q\ clS Qb function of time for different 
values of J. The field dynamics for J = 0.01 is already found to approximate 
well the limit J — ?■ 0"^, while for J = the classical field average vanishes as it 
should by symmetry. Similar to what is observed for the two-point function, 
0ci exhibits three characteristic time regimes. After an initial power-law 
growth, the dynamics slows down considerably. Finally, 0ci saturates to the 
order-parameter value in the spontaneously broken phase (see Fig. [3]). The 
field average takes on its maximum value when all configurations belong to 
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Figure 4: Upper graph: The equal-time two-point correlator as a function of 
time, spatially Fourier transformed and for different sources J with = 1, 
7 = 0.4. Lower graph: The classical field average as a function of time. For 
identically vanishing source (J = 0), 0ci remains zero by symmetry. 



the same phase and the phase transition is completed. No kink contributions, 
which interpolate between different field values, are expected in this case. 
Accordingly, by comparison with the upper graph one observes that the 
time when 0^1 reaches its maximum is also the time when the equal-time Fc\ 
exits the quasi- stationary period and starts to drop rapidly. A very similar 
behavior can also be observed for N = 2. 

It remains to show that the intermediate quasi-stationary period for the 
two-point function is indeed characterized by topological defects. Fig. \5\ 
shows the correlator as a function of momentum for fixed time and J = 0. 
The upper graphs are for 7 = 0.4 and t = 40. At this time the evolution of 
the two-point function is in the quasi-stationary regime both for A^ = 1 (left) 
and N = 2 (right). The fits using the kink ansatz (fT2|) with n = 0.06 
and (IT^ with ^ = 1.3 for textures are also shown. Simulation results 
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Figure 5: Classical simulation results for the equal-time two-point function 
and fits to the corresponding defect ansatze (solid lines) are shown as a 
function of momentum. Upper graphs are for 7 = 0.4 and lower graphs for 
7 = 0.02 both for = 1 (left) and N = 2 (right). The times are chosen to be 
within the quasi- stationary period of the two-point function. Data and fits 
agree well for lower momenta, indicating the presence of topological defects. 



and analytical expressions agree well in the infrared, clearly indicating the 
presence of topological defects. For comparison, we give in addition results 
for smaller damping, 7 = 0.02. In this case the evolution for the correlator 
enters the quasi- stationary regime for the equal-time two-point function at 
a later time beginning around t = 100. The lower graphs show the results 
at t = 1000, which again indicate the presence of defects. The classical- 
statistical simulation results will be compared with those from the 2PI 
effective action below, where we will also further discuss the J-dependence. 
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4 Quantum theory: 2PI to NLO 



To compute the quantum equations of motion for our damped model of 
Sec. 121 it is useful to note that the classical field equation (|2]) can also be 
viewed as the equation of motion of an undamped 2+1 dimensional scalar 
field theory in an expanding anisotropic geometry with the expansion factor 
a(t) = qq exp(7t) such that 'j = a/a This allows one to use standard 2PI 
effective action techniques for the quantum theory, which do not suffer from 
secularity problems and where nonperturbative approximations are available 
such as the 2PI large- expansion to NLO [5]. Further details can also be 
found in appendix |A1 

We consider the 0(iV)-symmetric quantum field theory with quartic self- 
interaction corresponding to the classical model of Sec. [2l For Heisenberg 
field operators ipa(t,x), we compute the quantum equations of motion for 
one- and two-point functions 



Here denotes the field expectation value, p is the spectral function 
determined by the commutator and F is the statistical propagator obtained 
from the anti-commutator of two field operators. Even though we are finally 
interested in spatially homogeneous systems, we consider for a moment the 
equations of motion for inhomogeneous fields and correlation functions to be 
used in Sec. O 

As detailed in appendix |Xl at NLO in the 2PI 1/A^ expansion one obtains 
the field evolution equation [5]: 



where the initial time is taken as t = and 0^ = 0a0a- Compared to the 
classical equation of motion ([2]), the effect of fluctuations in the quantum 
description enter in terms of the mass-like term and self-energy which 



= {ipa{t,x)), (16) 

Pabit,t';x,y) = i{[(pa{t,x),ipb{t',y)]) , (17) 
Fab{t,t']x,y) = -{{(pait,x),ipbit',y)}) - (j)a{t,x)(l)b{t',y) . (18) 




(19) 
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are given below. The corresponding evolution equations for F and p are [5] 

[{dl + ^dt - dl) 6ab + MUt, x)] FUt, t'- X, y) 

dt" [ dz a{t")T.P^,{t, t"; X, z)F,,{t" , t'; z, y) 



+ j dt" j dz a(r)Sf,(t, t"- X, z)pUt'\ t'- z, y) , 
[{d^^ + -idt - dl) 5ab + Ml,{t, x)] pUt, t'; X, y) 
= - I J t"; X, z)pUt'\ t'- z, y) . (20) 

At NLO in the 2PI expansion the mass and self-energy terms appearing 
in ([T9D and (EOD are [5] 



ab 



Ml,{t,x) = {p'' + -^[F,,{t,t-x,x) + <p\t,x)\y, 

+ ^ [Fab{t, t; X, X) + (j)a{t, X)(j)b{t, X)] , 

^ab{t^t']x,y) = -^(^lF{t,t']x,y)[(l)a{t,x)(l)b{t',y) + Fab{t,t';x,y)] 

1 

4 

1 

4 

Kbit^t'iX,y) = -^(^Ip{t,t';x,y)[4)a{t,x)(f)b{t',y) + Fab{t,t';x,y)] 
+lF{t, t'; X, y)pab{t, t'; X, y) + Pp{t, t'; x, y)Fab{t, t'; x, y) 
+PF{t,t';x,y)pabit,t';x,y)^ , (21) 



--/p(t, t'; X, i/)pa6(i, t'; X, y) + Pp(t, t'; x, y)Fab{t, t'; x, y) 
i'; a;, 2/)Pa6(t, t'] X, y) 
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with the summation functions 

A 



lF{t,t'; x,y) 



6N 



F\t,t'-x,y)-^p\t,t'-x,y) 



dt" I dz a{t")Ip{t,t";x,z) 



+2 / dt" / dz a{t")lF{t, t"; X, z)Fab{t", t'; z, y)pab{t", t'; z, y) 



Ip{t,t';x,y) = -^Fab{t,t';x,y)pab{t,t';x,y) 



X 



3N 



dt" / dz ait")Ip{t,t";x, z)Fab{t",t'; z,y)pab{t",t'; z,y) 



where F^ = FabFab, etc. For non- vanishing field one also has 

PF{t,t';x,y) = -^(^HFit,t';x,y) - dt" J dz a{t")[Hp{t,t"; x, z) 

xlF{t",t';z,y) + Ip{t,t";x,z)HF{t",t';z,y)] + j dt" jdza{t") 
X [HFit, t"; X, z)Ip{t", t'; z, y) + hit, t"; x, z)Hp{t", t'; z, y)] 

- I dt" I dt"' [ dzdu a{t")a{t"')Ip{t, t"- X, z)HF{t", t"'- z, v)Ip{t"\ t'- v, y) 
Jo Jo J 

+ [ dt" [ dt"' [ dzdu a{t")a{t"')Ip{t, t"; x, z)Hp{t", t'"; z, u)lF{t'", t'- z/, y) 
Jo Jo J 



dt" j^^ dt'" j dzdua{t")a{t"')lF{t, t"; x, z)Hp{t", t'"; z, u)Ip{t'", t'; u, y) 



Pp{t,t';x,y) 



3N 



Hp(t, t'; X, y) — dt" j dz a{t") 



X [Hp{t, t"; X, z)Ip{t", t'- z, y) + Jp(t, t"- x, z)Hp{t", t'; z, y)] 
+ [ dt" [ dt'" [ dzdu a{t")a{t"')I pit, t"; x, z)Hp{t", t'"; z, u)Ip{t'", t'; u, y) 
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with 



HF{t,t']x,y) = -(t)a{t,x)Fab{t,t']x,y)(l)b{t' ,y) , 
Hp{t,t';x,y) = -Mt,x)pab{t,t';x,y)Mt',y)- (22) 

We note that the terms —\lpP and —\PpP in (1^ and the terms proportional 
to in the summation function Ip describe genuine quantum corrections, 
whereas all other terms would also be present in a classical-statistical 
description of the dynamics jH]. Neglecting these terms one obtains, 
therefore, the evolution equations for the corresponding classical-statistical 
field theory [5]. 

The above evolution equations have to be supplemented by initial 
conditions for the field x) and its first derivatives, and similarly for the 
correlation functions. For the spectral function these are determined by the 
boson field commutation relations, i.e. 

pab{t,t']X,y)\t=t' = 0, 

Six — v) 

dtpab{t,t';X,y)\t=t' = 77^ Sab, (23) 

a{t) 

dtdt'pabit,t';x,y)\t=t' = 0, 

which are valid at all equal times including t = t' = 0. The initial conditions 
for the statistical propagator for homogeneous systems are chosen to be of 
the free-field form in accordance with ( fT5l) of Sec. E) 



Fab{t,t';p)\t ^ 



t=t' 


=0 


- „ Sab 

Zujp 


t=t' 


=0 


= 0, 








t=t' 


=0 


— —Oab, 







dtFab{t,t';p)\t=t'=o = 0, (24) 
dtdeFab{t,t';p)\t 

where again Up = \/p'^ + rn? with m? = —p^. The initial field average and 
also its initial time derivative are taken to vanish as in (fT4l). 



= 0) = , dt<Pa{t)\t=^ = 0. (25) 

By considering 0(A^)-symmetric initial conditions as above, the evolution 
equations can be simplified. Rotating the field appropriately we are only left 
with one non-vanishing field component and the propagators reduce to one 
parallel and (A^ — 1) times the same perpendicular component. 
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Figure 6: Comparison of the evolution of the field from the 2PI 1/A^ 
expansion to NLO (lines) and classical-statistical simulations (symbols) for 
= 1 and different values of the source J. 



For comparison of results from the quantum evolution equations and from 
the classical-statistical simulations of Sec. El we regularize the theory on a 
spatial lattice with periodic boundary conditions and consider discretized 
time steps At. Typically At = 0.1 is sufficient to observe insensitivity of 
results to time-discretization errors. We will chose the same parameters as 
in Sec. |3]if not stated otherwise. 



In the classical-statistical simulations of Sec. E] we considered different 



inhomogeneous and after averaging over many members of a given ensemble 
the homogeneous field average is obtained. In contrast, the quantum 
description of Sec. H] evolves the field expectation value using ( !T9|) . which 
includes already the quantum fluctuations and requires no further averaging. 
Therefore, for spatially homogeneous quantum systems the time-dependent 
field (f){t) has no spatial dependence. It is pointed out in Ref. [TB] that in 
the presence of topological defects an approximate description based on a 
1/A^ expansion can be inaccurate in this case. The results are based on a 
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homogeneous fields 



individual realizations of the classical field. 



Each realization is 
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comparison of classical-statistical simulations as in Sec. |3]with solutions of the 
NLO evolution equations fl2Up in the classical-statistical field theory limit as 
described in Sec. |H In Ref. [16] the field expectation value is taken to vanish. 
However, = does not correspond to the true field value in a symmetry 
breaking phase transition. In particular, the evolution equations fl20|) for 
the correlation functions F and p receive corrections from a non-zero cf). In 
this section, we extend that study by following the time evolution of the non- 
zero order parameter and for non-zero source J. After analyzing in some 
detail the shortcomings of the 2PI effective action approach for homogeneous 
fields to describe physics of topological defects, we proceed in Sec. [6] with an 
inhomogeneous 2PI description resolving these difficulties. 

We compute the nonequilibrium time evolution for the quantum system 
solving f lT9|) - fl22|) using homogeneity and isotropy, such that the field 
4>(t,x) = (j){t) and the statistical correlator F{t,t'; x,y) = F(t,t'] \x — y\), and 
equivalently for the spectral function p. The resulting equations are solved 
in spatial Fourier space, with the initial conditions (j){t = 0) = dt(f){t)\t=o = 
as well as f l2^ and for the correlators. We have explicitly verified from 
the numerical solutions that classical-statistical fluctuations dominate over 
quantum fluctuations for all times that are relevant for topological defect 
formationj^ In this case all terms dominate over terms in the above 
evolution equations. As mentioned in Sec. H] this represents the classical- 
statistical limit of these equations [H], which we employ in the following. 
As a consequence, the 2PI approach and the 'exact' classical-statistical 
simulations describe the very same theory and differences arise solely because 
of approximations. 

Fig. |6] shows the time evolution of the field (f){t) for = 1 and various 
values of the source J. The results from the 2PI 1/A^ expansion to NLO 
are drawn as lines while those from the corresponding classical-statistical 
simulation are represented by symbols. According to the discussion of Sec. [3l 
for our model we expect significant contributions from defects starting around 
t > 10 after which the equal-time two-point function F{t, t; p) enters a quasi- 
stationary regime. This quasi- stationary period for F{t,t;p) is absent for 
sources exceeding J* ~ 0.4 and lasts until about t = 40 for J = 0.04 (see 
Fig. H]). We plot in Fig. [6] the field evolution until t = 80. One observes 
that the initial power-law growth until t ^ 10 is accurately reproduced by 

^Of course, genuine quantum corrections are crucial for the asymptotic late-time 
approach to thermal equilibrium, which has been studied elsewhere [SJ IHl IS] • 
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Figure 7: The statistical two-point function as a function of momentum for 
fixed times t = 10 (top) and t = 40 (down) for = 1 and different values 
of the source J. Compared are NLO results (lines) and those from classical- 
statistical simulations (symbols). 
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Figure 8: Shown is the non-zero component of field expectation value as in 
Fig. E] but for N = 2. 



the NLO approximation. A rather good agreement is also found at later 
times for J = 0.4, for which defects are expected to play no role. This 
agreement is remarkable by itself since there is no a priori reason that a 
1/N expansion evaluated at = 1 should be accurate. In contrast, for 
smaller J significant deviations between classical-statistical simulation and 
NLO results are building up quickly. 

A similar finding can also be obtained from comparisons of the results 
for the statistical two-point function. Fig. [7] shows F{t,t]p) as a function of 
spatial momentum for fixed times and various values of J for A^ = 1. The 
upper graph at t = 10 confirms a rather good agreement between the NLO 
approximation and simulation results for small enough J. Only for J = 0.4 
some deviations are observed at t = 10 at which F{t, t; p) started already 
to decay significantly according to Fig. |4] of Sec. |3l In this region F{t, t; p) 
oscillates rapidly, so a slight phase difference between the two computations 
can lead to a pronounced difference. At later times this dynamics slows 
down and the observed deviations become smaller again. This can be seen 
from the lower graph of Fig. [7] showing results at t = 40. In contrast, for 
smaller J one observes increased deviations at t = 40 as compared to earlier 
times. Comparison with Fig. Oof Sec. [3] shows that the NLO results do not 
reproduce well the shape of F{t, t]p) which is characteristic for contributions 
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Figure 9: The same as in Fig. [7] for N = 2. Shown is the longitudinal part 
of the statistical two-point function. 
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Figure 10: The same as in Fig. [9] but for the transverse part of the statistical 
two-point function. 
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from defects for small J. 

The same analysis is also done for N = 2. The qualitatively similar 
results are shown in Figs. |8] to [TOl In general, we find a somewhat 
better agreement between the NLO approximation and classical-statistical 
simulations as compared to = 1, which should be expected for a 1/N 
expansion. We have verified that for N > 2 the agreement becomes excellent, 
as is reported as well in Ref. [TB]. However, no defects are present in this case. 
In view of these results, it would be very interesting to perform a comparison 
with the 2PI expansion to NNLO. However, the computational effort 
becomes considerable at NNLO [15] and alternative approaches may be more 
appropriate, which is discussed in the following. 

6 2PI expansions with inhomogeneous fields 

In Sec. [5] the 2PI effective action is considered for a time-dependent but 
spatially homogeneous field expectation value This reflects the spatial 
homogeneity of the considered quantum system. Apparently, this is an 
inefficient starting point for a 1/A^ expansion if physics of inhomogeneous 
classical field configurations such as topological defects is involved. In 
classical-statistical simulations each realization of the field is inhomogeneous 
and only after averaging over many members of a given ensemble the 
homogeneous field is obtained. In contrast, the quantum 2PI effective 
action does not give information about individual realizations but only about 
expectation values or quantum-statistical averages. It is also not obvious 
that one could consider individual evolutions of inhomogeneous quantum 
field expectation values whose ensemble average would give the correct 
homogeneous field without " double counting" . 

Here we propose an approach, which takes into account the fact 
that defect formation is most efficiently described using inhomogeneous 
fields but without the need of sampling techniques such that no double 
counting problem for quantum theories can occur. In this prescription the 
spatially homogeneous quantum system is constructed from a corresponding 
inhomogeneous system and a volume average. More precisely, the 
homogeneous field is replaced by an inhomogeneous field x) with 
an additional spatial x-dependence. Similarly, the spatially homogeneous 
two-point correlation functions such as F{t, t'; x — y) are represented by their 
inhomogeneous counterparts, F(t,t'; x,y), etc. The inhomogeneous field and 
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correlation functions are governed by the quantum evolution equations fll9p - 
(122]) for the 2PI expansion to NLO. The original homogeneous field 
and correlation functions are obtained from averaging the corresponding 
inhomogeneous ones over some local volume V. In general, this "coarse 
graining" volume V is smaller than the system size L'^ and is determined 
by the statistical error one wants to achieve. It can also be convenient to 
implement the averaging in Fourier space, which we employ below. 

In the following we show that the inhomogeneous 2PI approach contains 
the classical-statistical simulation results of Sec. [3] already at leading order 
(LO) in the large- expansion. The equations are given by (fT9|) - (l2T|) with 
S-^ = = at LO and read 

(^d^ +jdt-dl+fi^ + -^ [F,,{t, t- X, X) + <P\t, X)] ^ Ct>a{t, X) = Ja{t, x) , 

(^d^ + ^dt-dl + fi' + -^ [F,,{t, t; X, x) + (P\t, x)] ^ Fab{t, t'; x,y) = 0, 

dl +^dt~dl + ii^ + -^ [F,,{t, t- X, x) + 02 (t, x)] ^ pab{t, t'- X, y) = . 

(26) 

One observes that the dynamics for the spectral function p does not influence 
the evolution of the field and the statistical function F at this order. 
We will not consider it in the following. Related dynamics based on an 
inhomogeneous Hartree approximation have been considered in Ref. [T7] . 

The inhomogeneous evolution equations fl26l) have to be supplemented by 
initial conditions, which we formulate here in spatial Fourier space using the 
functions 

Xa(t = 0,i.) = -i=e^-(^) , dtXa{t.p)\t=o = ^e^^^^\ (27) 

again all in units of appropriate powers of a^, with random numbers aa{p) 
and Pa{p)- From this we compute the Fourier transform to obtain Xa{t = 0, x) 
and derivatives. The initial correlation functions are given by (!23|) for the 
spectral function and for the statistical function we take 

Fabit,t';x,y)t=t'=o = Xa(0, x)x6(0, y) , 
dtFab{t,t';x,y)t=t'=o = 0, (28) 
dtdt>Fab{t,t';x,y)t=t'=o = dtXa{t, x)t=odt'Xb{t' ,y)t'=o , 
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from which we also determine the initial spatial derivatives. The initial field 
and also its initial derivatives are taken to vanish with 



This corresponds to the Gaussian initial conditions for averaged quantities 
as employed in f lT^ and f lT5|l . or fl23l) - fl25l) of previous sections. 

We consider J = for which the field vanishes at all times, (f)a{t,x) = 0, 
with the zero initial field fl29|) and derivatives. Using fl28|) one directly verifies 
the equivalence of the LO evolution equation for F with the classical field 
equation ([2]) if the latter is multiplied by one power of the field. As a 
consequence, at LO the dynamics is purely classical. It should be emphasized 
that the initial conditions are the same only on average. Therefore, the 
differences concern the inhomogeneous initial conditions and the averaging 
procedure to obtain the homogeneous correlation functions: While in the 
classical-statistical simulations of Sec. [3] the average is done over many 
individual runs with inhomogeneous fields, here the result is obtained from 
a single integration of the differential equation for Fabit,t'] x,y) and a 
volume average. In this case, the comparison between the classical-statistical 
dynamics of Sec. |3]and the present LO dynamics boils down to a comparison 
between ensemble and (local) volume averages. 

To demonstrate that the volume average of F{t,t'] x,y) indeed becomes 
homogeneous and that the results agree well to those obtained from the 
classical-statistical simulations of Sec. |3l we integrate the evolution equation 
for F{t, t'; X, y) for the = 1 component theory in configuration space on 
a spatial lattice with L = 8192 and periodic boundary conditions. We 
then spatially Fourier transform the results to get correlation functions 
in momentum space, in which we also do the averaging to obtain 
the homogeneous correlation function. More precisely, we evaluate the 
inhomogeneous two-point function F{t,t']p,p') at p' = —p and compute the 
homogeneous average from binning F{t,t']p,—p) for all p belonging to a 
given bin on a logarithmic scale. In Fig. [TT] we show the result using 1024 
bins, along with the statistical error. We plot the correlation function as a 
function of momentum at t = 40. At this time important contributions from 
topological defects can be observed according to the discussion of Sec. |3l For 
comparison we also give the previous classical-statistical simulation result 
and the fit for the analytical defect ansatz. The agreement is very good and 
the statistical errors barely visible from the plot. This demonstrates that 




(29) 
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Figure 11: Comparison of results for the statistical two-point function 
as a function of the momentum at t = 40. The results from the 
inhomogeneous 2PI large-iV expansion at LO (symbols) and from classical- 
statistical simulations (dotted line) agree very well. Shown is also the kink 
fit (solid line) of Sec. [3] 



the inhomogeneous 2PI approach contains the classical-statistical simulation 
results of Sec. [3] already at LO in the large- expansion. 

At NLO in the 2PI 1/A^ expansion genuine quantum corrections enter 
the dynamics. Here it is important to note that in this approach, classical- 
statistical fluctuations and dynamics is included exactly already at LO, 
whereas quantum corrections are treated in an approximate, truncated form. 
Though the direct integration of the corresponding NLO equations (|T9ll ~ 
(l22l) along the same lines is straightforward in principle, it is rather costly 
numerically, in particular in higher dimensions such as d = 3. We defer 
the discussion of suitable simplifications, which allow one to describe defect 
formation in phenomenologically relevant quantum field theories in practice, 
to a subsequent publication [20] . 



Much of our theoretical knowledge about the dynamics in complex physical 
situations such as defect formation, nonequilibrium instabilities or turbulence 



7 Conclusion 



24 



relies on classical simulations. Classical-statistical simulations are based 
on inhomogeneous field dynamics and importance sampling to describe 
homogeneous systems. In many cases these are known to be reliable 
approximations for quantum systems, at least for short times before the 
approach to thermal equilibrium characterized by Bose-Einstein or Fermi- 
Dirac distributions. 

In this work we have proposed a quantum approach which combines the 
successful aspects of inhomogeneous classical simulations with the ability to 
take into account quantum corrections. This is necessary since standard 
nonperturbative quantum approaches, such as the 2PI large-A^ expansion 
for homogeneous fields, turn out to be not appropriate for a quantitative 
description of defects for one- and two-component theories at low expansion 
orders. We have demonstrated that the result of Ref. [16] is correct, and 
that the discrepancy is indeed qualitative rather than quantitative. To see 
this we considered the procedure of introducing a source to the equation and 
continuously bringing the source to zero. At some threshold value of J, the 
potential acquires degenerate minima, thus allowing for defects, and this is 
the point where the discrepancy is seen to kick in. Introducing a non-zero J 
also leads to including a non-zero field expectation value, which is necessary 
to follow the true field value in a symmetry breaking phase transition. The 
omission of this in Ref. [16] was one of the main motivations to start this 
investigation. In general, we find that NLO results for N = 2 agree somewhat 
better to those obtained from classical-statistical simulations than for = 1 
as expected for a large- iV expansion. It would be very interesting to extend 
the study to NNLO, however, a more appropriate calculation should be 
based on the 2PI effective action for inhomogeneous fields and a volume 
average. We emphasize that the latter does not involve any double counting 
which could appear in sampling prescriptions for inhomogeneous quantum 
evolutions. We have shown that already at leading order in a 1 /TV expansion 
the results from classical-statistical simulations are recovered. At NLO 
quantum corrections enter the dynamics. This provides a clear procedure to 
describe defect formation in quantum field theory, though it is numerically 
very costly. The computations of the inhomogeneous NLO approximation 
along with a discussion of suitable simplifications will be presented in a 
subsequent publication [20] . 

A further promising application of our method concerns the question of 
thermalization in quantum chromo dynamics in the context of heavy-ion col- 
lisions. Classical-statistical simulations in lattice gauge theory have been 
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applied to plasma instabilities and the question of isotropization at early 
times 122] • At low momenta the instability dynamics is expected to be 
well approximated by classical methods. However, at higher momenta quan- 
tum corrections may play an important role, because occupation numbers are 
low at large momenta. A 2PI loop expansion for inhomogeneous gauge fields 
along the lines of this work would include the classical-statistical simulation 
results already at lowest order. Going beyond lowest order in a 2PI loop ex- 
pansion may be an appropriate way to include relevant quantum corrections 
at higher momenta. 

This work is supported in part by the BMBF grant 06DA9018. 

A Calculation of the 2PI effective action 

The classical field equation (|2]) contains a damping term, which can be 
obtained from a principle of least action if it is viewed as the equation of 
motion for an undamped 2+1 dimensional scalar field theory in an expanding 
anisotropic geometry with metric [H] 

ds'^ = -df + dx^ + a^{t)dy^ , a{t) = cq exp(7t) . (30) 

Then (j2]) follows from the equation of motion for the field (j){t,x,y), 

dUa + -dt(l>a - dl<j)a ' ^d^a + /i^a + -^AMh)<Pa - Ja = (3l) 

by neglecting the dependence of the field on the spatial variable y. 

A.l 2PI effective action in expanding geometry 

To obtain a quantum description of our model, we will provide some details of 
the calculation of the 2PI effective action in an expanding geometry in this 
section. For similar discussions see Refs. [TJ [23]. To keep the formulation 
general, we consider an arbitrary metric g^'^ with signature ( — h + . . . ). For 
a compact notation variables as x = (a;°,x) are taken to denote space-time 
variables in c? + 1 dimension, in contrast to the convention in the rest of this 
paper. At some points we insert the special metric of our model to recover 
the equations of Sec. HI For a general metric g^'^ the classical action reads 

S[v; 9]= [ ^/^) (-]-g^''{x)d,Va{x)d.^a{x) - V{^') \ . (32) 
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Here = J,^ dx^ J d'^x involves integration over the Schwinger-Keldysh 
closed time contour [T^ . The determinant of the metric tensor is denoted 
as f7 = detg^'^. For our model discussed in this work we have 

^/—g{x) = a(x°) = ao exp (70;°) . (33) 

The inverse classical propagator is given by 



-9{x) ^/-giy) S<^a{x)6^b{y) 
-Dab{x;if) — , (34) 

V-9{y) 



where for our model D reads 



(pa{x)(pbix). (35) 



3N 

Introducing sources Ja{x) and Rab{x,y) we construct the generating 
functional 

Z[J,R;g] = j ^ipexp!^i(^S[(p]g] + j y^^g{x)Ja{x)(pa{x) 

+\ j V-9{x) \/-g{y)Rab{.x, y)^a{x)^b{y)^ | , (36) 
W[J,R;g] = -i\nZ[J,R-g], (37) 

where W[J, R; g] is the generating functional for connected Green's functions. 
The field expectation value (j)a and the connected two-point function Gab are 
defined as 



^/-g{x) ^Ja{x) 
1 1 5W[J,R]g] _ 1 

^-g{x) ^-g{y) ^Rab{x,y) 2 



Ux)<Pb{y) + Gab{x,y)] . (38) 
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The 2PI effective action T is the Legendre transform, 

run ^ um ^ f ^jniR^ f SW[J,R;g] ^ , 
r ^,^;^ = W\J,R\g\ - Ja{x) - / — — -Rab{x,y) 

= W[J,R;g]- I ^-g{x)ct>a{x)Ja{x)- \ ! ^-g{x)^-g{y) 

J X ^ J x,y 

xRab{x, y) [(t)a{x)(t)b{y) + Gab{x, y)] , (39) 
which can be written as [18] 



r[0, G- g] = S[<P- g] + ^Tr In G^^ + '-l^iG^^G + V^i^P, G; g] + const. , (40) 

where r2 contains all two-particle irreducible contributions. The equations 
of motion for the field expectation value and the propagator are 

y ^ , , = -Ja{x)- [ V-9{y)Rab{x,y)^b{y) , 

^J-g{x) 0(l)a[x) Jy 

^J-g{x) ^-g{y) SGab{x,y) 2 

For Rab{x,y) = the second equation corresponds to the Schwinger-Dyson 
equation 

G-b\x, y) = G^l^ix, y- ^) - S,,(a:, y) , (42) 
where the self-energy E is defined by 

1 1 6T,[4>,G;g] 



^ab{x,y) = 2i 



V-9ix) V-9iy) ^Gab{c,y) 
<o) - y) , - 



= -^El^^^(a:)^^=^ + S„,(x,y) (43) 

In the last equation we have separated a local part S^°''(x) from the self- 
energy, which for our model is given by 

Sl?(^) = ^F,c{x,x)5,b + ^^a.(x,x) . (44) 



Multiplication of (1421) with G and integration leads to the differential 
evolution equation for the propagator 

D^,(a;;0) + Sl°)(a;)l Gb.{x,y) 



: / ^/-g{z)T.ab{x,z)Gbc{z,y) - iSac ^ . ^ ■ (45) 

v-9{y) 
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A.2 2PI l/N expansion to NLO 

For the 0(A^)-symmetric scalar field theory with classical potential (j3]) we 
classify the contributions to the 2PI action by their scaling in powers of 

[0, G- g] = r^o [G; g] + o [<f>,G;g] + .... (46) 

The leading order contributions (LO) scale as A^^, the next-to- leading order 
terms (NLO) by and so on. According to Ref. [S] at LO this is given by 

r^O[G; g] = ^ ^/^^Gaaix, x)G,,(x, x) (47) 

and the NLO contributions are 



rfO[0,G;^] = ^Trln[5(G)] 
^ 6N 



ix)\/-g{y)Iix, y)(t)a{.x)Gab{,x, y)(f)biy) (48) 



with the definition 



B{x, y) = Jl^ + ^^^Gaf>{x, y)Gat{x, y) (49) 



and the integral equation 



H^^y) = -^Gabix,y)Gabix,y) 



- \/-9{z)I{x, z)Gab{z, y)Gab{z, y) . (50) 

In order to write the evolution equation (145|1 in a form which is suitable 
for numerical computations, we identically decompose the time-ordered 
propagator in a statistical part F and a spectral part p [3], 



Gab{x, y) = Fabix, y) - -sgn,^(a;° - y^)pab{x, y) . (51) 
Similarly, for the non-local part of the self-energy defined in (H3i) we employ 

Safe(x, y) = ^ab{x, y) - -sgn.^(x° - y°)S^^(x, y). (52) 
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Applying this decomposition to the evolution equation ( I45p and separating 
the real and the imaginary part leads to the evolution equations for F and p 



Fhc{.x, y) 

dz' j d''z./^^El,{x,z)F,,{z,y) 
+ dz^ j d'^z^y-g{z)T,^h{x, z)pbc{z, y) , 



Pbc{x,y) 

dz° I d'^z^/-g{z)T.^^f^{x, z)pbc{z, y) . 



(53) 



From equation (14T1) with vanishing quadratic source R we can derive the 
corresponding evolution equation for the field, 



A 



- Ja{x) 

dy' I dS^^MKb{^^yA = mb{y)- 



(54) 



A. 3 Energy- momentum tensor 



The energy-momentum tensor can be calculated by varying the 2PI effective 
action with respect to g^'^{x), 



2 ST[4>,G-g] 



(55) 
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At NLO in the 2PI expansion this gives 

dXFaa{x,x') - \g,, [g^^dldtda, + M, 



\9^''dMx)dx(t)a{x) + ]^^?(t)a{x)(t)a{x) + ^ (0, (x)</), (x) ) ' 



+ 



Fab{x,x') 



A 2 1 

-9f^y{x)^ [Faa{x,x)] - g^^{x)-lF{x,x) 



-9f,u{x) 



Pf{x,x) + :^HFix,x) 



(56) 



with the summation functions If, Pf and Hf as defined in Sec. HI 
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